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Each and every energy dependent partial-wave analysis is parameterizing the pole positions in 
a procedure defined by the way how the continuous energy dependence is implemented. These 
pole positions are, henceforth, inherently model dependent. To reduce this model dependence, we 
use only one, coupled-channel, unitary, fully analytic method based on the isobar approximation 
to extract the pole positions from the each available member of the worldwide collection of partial 
£N) ' wave amplitudes which are understood as nothing more but a good energy dependent representation 

of genuine experimental numbers assembled in a form of partial-wave data. In that way, the model 
dependence related to the different assumptions on the analytic form of the partial-wave amplitudes 
is avoided, and the true confidence limit for the existence of a particular resonant state, at least in 
one model, is established. The way how the method works, and first results are demonstrated for 
the Sn partial wave. 

PACS numbers: 14.20.Gk, 12.38.-t, 13.75.-n, 25.80.Ek, 13.85.Fb, 14.40.Aq 
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When resonances are associated with the eigenstates of the complete Hamiltonian for which there are 
only asymptotically outgoing waves, their identification with scattering theory poles is unquestionable. This 
statement is in details elucidated in ref. Consequently, in order to get the full information about physical 
systems and resonant states under observation, we must be entirely focused onto analyzing and interpreting 
the scattering matrix singularities of the Mandelstam analytic function 0] obtained from experiments. While 
, the value of the scattering amplitude on the positive energy cut defines the physical amplitude in the s or u 

channel depending whether we approach the physical axes from above or below, the simple poles which are 
^vq | situated on the physical axes in a subthreshold region are related to the bound states. As it is believed that 

there is no fundamental difference between a bound state and a resonance, other than the matter of stability, 
when simple poles of the coupled channel amplitude occur on unphysical sheets in the complex energy plane, 
they are to be associated with resonant states [|| . 
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INTRODUCTION 



The fact that we are trying to extract the value of a quantity lying in the complex energy plane while 
performing experiments only on the physical axes, is the essence of all problems, and origin of many mis- 
\ understandings occurring in the literature. Namely, each pole is not only squatting in an experimentally 

inapproachable domain, but is simultaneously governing each and every process between all allowed few body 
channels. However, we usually measure observables only in one channel at a time. If the single-channel 
observables are measured, we obtain the single- channel scattering amplitude, and we only get the pole positions 
in one channel. Nevertheless, due to the Mandelstam hypothesis, these poles are affecting all channels, so we 
have to treat them all and not just the measured one. Consequently, the underlying theory which is to be used 
to find the scattering matrix amplitude must be a coupled-channel one, and of course analytic and unitary. 
And this is not the end. Once we have found the coupled-channel scattering matrix amplitude, we have to find 
and quantify all its poles. Unfortunately, this is not a simple task, so each partial wave analyses, even in a 
multi-channel case, has its own way of parameterizing this inaccessible quantity. The result is that the model 
dependencies are introduced. 



This brings us to various ways how the complex energy plane poles are up to the present moment parame- 
terized in the literature. First attempts are done with single-channel partial wave amplitudes, and the oldest 
and most frequently met way is the concept of Breit-Wigner parameters. 
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The initial attempts to use the Brcit-Wigncr function with constant parameters to represent the scattering 
matrix amplitudes on the physical axes immediately revealed the fact that this function is too simple. More 
terms were needed. One had to introduce the energy dependent background, and one had to do it in a unitary 
way. Unfortunately, for quite some time it has been known, but not commonly accepted, that a unitary addition 
of background terms influences the peak position of the scattering matrix absolute value on the real axes in spite 
of the fact that the pole position is not changed. Peak position is an interplay of Breit-Wigner parameters and 
background terms. And the peak position is the quantity which is usually extracted from experiments. Conse- 
quently, when Brcit-Wigncr parameters defined in such a manner are chosen to represent the pole position, they 
must be background dependent, and the only case when the Breit-Wigner parameters do exactly correspond to 
the pole position is when we have accidently guessed the correct form of the energy dependent background. If 
the background is wrong, Breit-Wigner parameters are not the pole parameters, but something else. And that 
is the reason why Breit-Wigner terms in general are not the pole positions, and are inherently model dependent. 

There are basically two ways to account for the background contributions. The first one, to unitary add 
energy dependent background terms to the constant-parameter Breit-Wigner function, is described afore. The 
second one is to allow the Breit-Wigner parameters to become energy dependent. That is predominantly done 
by modeling the Breit-Wigner width [4|t9|]- 

There is a number of ways to introduce energy dependent Breit-Wigner width. In ref. Q energy dependent 
width is a part of resonant term of theoretical function which is associated with the T-matrix near the 
resonance. In refs. [H, @ energy dependent width is related to the resonant part of the S-matrix. In a method 
proposed in ref Q and applied in refs. [H, Q width is dcffincd from the function consisting of a background 
term and Brcit-Wigncr shape term. 

One well known method for treating the nearby channel openings is the Flatte formula. The Flatte's method, 
introduced in 1976 [Toj . is recognizing the fact that the partial- wave T-matrix feels the presence of new channel 
openings, and it is taking it into account effectively. Flatte proposes to modify the traditional Breit-Wigner 
form by assuming that the width with becomes proportional to the phase space. The amplitude poles are then 
again represented as the singularities of the modified Brcit-Wigncr function. 

The fact that the Breit-Wigner terms in general are not the pole positions, and are inherently model 
dependent, was timidly mentioned by several authors (see for instance ref. [111]). That was first strongly 
pointed out by Hohler in refs. [lU, [l3j], where the definition of "local Breit-Wigner fit" and the concept of 
"searching for the pole position" using speed plot technique were introduced. Hohler clearly distinguished 
between Breit-Wigner parameters (which should in the absence of a better way be obtained by locally fitting 
partial wave amplitudes with a Breit-Wigner function plus some background terms) and pole parameters 
which should be obtained, as he recommended , by the speed plot technique. He has always been pointing out 
that Breit-Wigner parameters are model dependent, and he continuously objected to compare them directly. 
His last warning was published not so long ago [l4j . However, due to unclear historical reasons, the practice 
of direct comparing Brcit-Wigncr parameters coming from different origin continued in Particle Data Group 
(PDG) compilations. Breit-Wigner parameters, extracted with different background parameterizations are still 
directly compared [l5| , averages are made and error analysis is performed neglecting the fact that they may be 
in fact completely differently defined parameters. This practice should be abolished. 

There is a long history of efforts to avoid the concept of Brcit-Wigncr parameters, and to look directly for 
the genuine pole positions. 

The first, and most frequently met method, is the speed-plot technique introduced by Hohler for the 
single channel scattering amplitudes. It is based on the idea already mentioned in ref. Q that the pole 
position should be found by expanding the scattering amplitude in the vicinity of the pole, and the speed-plot 
technique is recommending to retain the first term only. This method is in principle acceptable if we are 
dealing with isolated poles far away from any nearby thresholds, but my fail otherwise. There is a number 
of cases where the methods can not be applied at all, and the best example was inability to use it to obtain 
the well known Sn(1535) resonance. The limitations of the method have been discussed by Ceci et al. [l6| 
where it has been shown that speed-plot technique is only the N=l term of a more general but demanding 
"regularization" method based on finding the N-th derivative of the scattering amplitude, and using it in a 
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local, three- parameter fit to the partial wave data (l7| . 

In the early fifties the time delay technique is introduced into scattering theory by several authors [lH-HH , in 
a way that they obtained expression for the time delay in a collision. Time delay, or in another words the time 
lapse between asymptotic states, can be directly related with phase shift of the T matrix. For further details 
on interrelation between speed plot and time delay see ref. [22j . 

The N/D method is a technique in which the dispersion relations are used to construct the amplitudes in the 
physical region using the knowledge of the left-hand cut singularities. The idea is to represent the partial-wave 
amplitude as a ratio of two functions. The numerator is represented with a function N(s) which is analytic in 
the s-plane only on the left-hand cut, and the function D(s) that is analytic on the right-hand cut only. The 
poles of the scattering amplitude are identified with the zeroes of the D(s), and the problem of extra zeroes is 
often difficult to be solved. The method has been introduced a long time ago by [23[, and since then it has been 
mostly used in meson physics, typically for cases when the knowledge about the left-hand cut is available [24l.l25l|. 

All enlisted methods arc good for the pole search within certain approximations, but yet we have to point out 
that the proper procedure to look for the scattering matrix poles is the full analytic continuation of scattering 
matrix amplitudes into the complex energy plane within a given model. 

In coupled-channel calculations the importance of the pole-search has recently been fully recognized. Some 
groups have offered more or less detailed concepts of their analytic continuation procedures [26|, [27j, while 
others have reported that the complexity of the analytic continuation of all Fcynman amplitudes of their model 
is beyond their reach [28}. Therefore, they had to rely on speed-plot technique entirely. In most cases, the 
analytic continuation procedure is rather cumbersome. 

The VPI/GWU collaboration clearly distinguishes the difference between Breit-Wigner parameters and 
pole positions, and states that "Poles and zeros have been found by continuing into the complex energy 
plane". Unfortunately, they fail to provide any details of their procedure. The EBAC collaboration also 
makes an analytic extrapolation of their amplitudes, and has recently presented a more detailed elaboration 
of their procedures [27l |. Other groups have extracted their pole positions using single-channel techniques 
such as speed-plot and time delay [28l - t3l| . Recognizing the importance of a direct analytic extrapolation, 



Dubna-Mainz-Taipci collaboration has recently performed the full analytic continuation, and in ref. 17| offered 
the reliable pole positions of their model. 



In spite of all these efforts, the question of systematic uncertainties still remains unanswered, because each 
model has its own, particular analytic form. So we wonder how stable, with respect to the analytic continuation 
procedure, these pole positions actually are. 

In order to get a reliable answer to this questions, we have decided to use only one method to extract pole 
positions from all published partial waves analyses, and compare the results. And we have chosen the T-matrix 
Carnagie-Melon-Berkeley (CMB) model. In other words, we take all sets of partial-wave amplitudes, treat 
them as nothing else but a good, energy dependent representations of all analyzed experimental data, and 
extract the poles which are required by the CMB method. In this manner, all errors due to different analytic 
continuations of different models arc avoided, and the only remaining error is the precision of CMB method 
itself. We shall also compare the obtained poles with the poles of each individual publication, and draw certain 
conclusions about features of individual methods as well. 



The general idea of this article is to recommend the possibility how to, in a maximally model independent 
way, simultaneously find all scattering matrix poles from the world-wide collection of partial wave amplitudes. 
We present the way of eliminating most systematic errors in analytic extrapolation by using only one, well 
defined procedure to extract pole positions for published partial wave amplitudes, understanding them as 
nothing more but a very confident energy dependent representation of all experimental data. 

To avoid congesting the reader with unnecessary information, wc in this paper in details illustrate how 
this method works for the Sn partial wave only. We show that N(1535) and N(1650) Sn resonant states are 
unambiguously seen in all analyzed PWA data, while the performed pole-search procedure strongly suggest 



4 



the existence of at least one more pole position in the vicinity of 1800 MeV. Therefore, all published PWA are 
consistent with the new Sn(1846) state seen in photo-production channel [2£|[32|]. We demonstrate that the 
existence of fourth Sn state around 2100 Mev is not excluded by any PWA, and is actually favored for the 
hadronic Dubna-Mainz- Taipei amplitudes [29|, |30[ • We compare the obtained results with the results published 
in literature, and make a final conclusion on the actual position of partial wave poles. 

However, the issue also arises how strongly the recommended pole-extraction procedure (CMB model) depends 
upon its own the model assumptions. Namely, CMB model has a number of assumptions, and it is very important 
to know how stable the pole positions are if CMB model choices are strongly modified. We have tested this 
problem extensively, and for the answer to this question we refer the reader to a parallel publication submitted 
to this journal (33j . 



II. FORMALISM 



The CMB model is isobar, coupled-channel, analytic, and unitary model, where the T matrix in a given 
channel is assumed to be a sum over the contributions from a number of intermediate particles (resonance 
and background contributions). The coupling of the channel asymptotic states to these intermediate particles 
determines the imaginary part of the channel function, and is represented effectively with a separable function. 
The real part of the channel function is calculated by the dispersion relation technique, thus ensuring analyticity. 
Besides the known resonance contributions, the background contributions are included via additional terms with 
poles below the 7rN threshold. Due to the clear analytic and separable structure of the model, finding the pole 
positions in CMB model is trimmed down to the generalization of the dispersion integral for the channel 
propagator from real axes to the full complex energy plane, and this is a very well defined procedure. In 
practice, we instead use a very stable, and numerically much faster analytic continuation method based on the 
Pietarienen expansion [42j in order to extrapolate the real valued channel propagator into the complex energy 
plane. 



A. Formulae 

Our current partial-wave analysis [34[ is based on the manifestly unitary, multichannel CMB approach 
of ref. [35l | . The most prominent property of this approach is analyticity of partial waves with respect to 
Mandelstam s variable. In every discussion of partial- wave poles, analyticity plays a crucial role since the poles 
are situated in a complex plane, away from physical region, and our measuring abilities are restricted to the 
real energy axis only. To gain any knowledge about the nature of partial- wave singularities would be impossible 
if partial waves were not analytic. Therefore, the ability to calculate pole positions is not just a benefit of the 
CMB model's analyticity but also a necessity for resonance extraction. In this approach, the resonance itself is 
considered to exist if there is an associated partial-wave pole in the "unphysical" sheet. 

We use the multichannel T matrix related to the scattering matrix S as: 

S a b(s) = 5 ab + 2iT ab (s), (1) 
where 5 a b is Kronecker delta symbol. The T-matrix matrix element is in the CMB model given as: 

where a(b) represents the outgoing (incoming) channel. In our analyses we use a,b = irN, r]N, ir 2 N. The 
initial and final channel b(a) couple through intermediate particles labeled i and j. The factors ji a are energy- 
independent parameters occurring graphically at the vertex between channel a and intermediate particle i and 
are determined by fitting procedure. Also occurring at each initial or final vertex is form factor fa L {s) 

tf ' ( "-( <w<a.+* )'' (2) 
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and phase-space factor p a (s) 

Pa{s) = (3) 

where s = W 2 is a Mandelstam variable, and q a (s) is the meson momentum for any of the three channels given 
as 



/ \ _ \/(f -(m + m a ) 2 )(s - (to - m a ) 2 ) 
9a(Sj- ^ • (4j 

Furthermore, L is the angular momentum in channel a. and Qi a , Q2a are constants. The factor f^ L (s) 
provides appropriate threshold behavior on the right-hand cut, and also produces a left-hand branch cut in the 
s plane. Parameters Q\ a and Qi a are chosen to determine the branch point and strength of the left-hand branch 
cut. In our analyses they have been taken to be the same, and are fixed to the mass of the channel meson a. 

Gjj L is a dressed propagator for partial wave JL and particles i and j, and may be written in terms of a 
diagonal bare propagator G\^ L and a self-energy matrix Ej^ using Dyson equation 

n jl 

G{ 3 L {s) = G^ L (s)+J2 G^ L ( S )Z{t( S )G^( S ). (5) 
k,i=i 



The bare propagator 



g ojl {s) = ^ (6) 
J s 



has a pole at the real value s^. The sign = ±1 must be chosen to be positive for poles above the elastic 
threshold which correspond to resonance. 

The nonresonant background is described by a meromorphic function, in most of the cases consisting of two 
terms of the form ([6]) with pole positions below it N threshold. For that case, the signs of the terms are opposite. 
The positive sign correspond to the repulsive and the negative sign to the attractive potential. In principle the 
number of poles can be increased arbitrarily (see the next subsection on background representation), but in 
reality the number is never larger than three. 

is the self-energy term for the particle propagator 

Ztt(s)=Y,rL L -<S> J a L (s)-lta L (7) 

a 

The (s) are called "channel propagators". They are constructed in an approximation that threats each 
channel as containing only two particles. We require that T J ^ have, in all channels, correct unitarity and 
analyticity properties consistent with a quasi- two- body approximation. 

The imaginary part of $^ L (s) is the effective phase-space factor for channel a: 

l^^i L {s) = [f J a L {s)fp a {s). (8) 

The channel propagator is evaluated on the real axes only 

r / \i2L + l 

***(»)= - [q[ X)] 0) 

2 , r„^M 2 



where by x we stress the fact that values are on the real axes. The real part of $^ L (x) is calculated using a 
subtracted dispersion relation 



(10) 
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FIG. 1: Parameterization of channel-intermediate particle vertex in CMB model. 



where x a = (m + m a ) 2 . For better understanding, the structure of the channel-intermediate particle form 
factor is given in Fig. [1] 

We give a matrix form of the final T matrix as defined in Eq. fl5J}: 



f(s) = \J 7m$(s) • 7 T • 



G (s) 



7-$( S )-fr .G (s) 



7 • y Jm$(s) 



(11) 



B. General idea 



In this paper we propose to use a method based on coupled-channel formalism, apply it to all partial wave 
data and partial wave amplitudes available "on the market" , and simultaneously analyze the underlying analytic 
structure. We have decided to use only one model to extract pole positions from all published partial waves 
analyses in order to evade the model assumptions of each approach, and compare the results on the same footing. 
And we have chosen the T-matrix Carnagic-Melon-Bcrkclcy (CMB) model. In other words, we take all sets 
of partial-wave amplitudes, accept them as nothing else but good representations of all analyzed experimental 
data, and extract the poles which are required by the CMB method. In this manner, all errors due to different 
analytic continuations of different models are avoided, and the only remaining error is the precision of CMB 
method itself (see ref. (HI). Of course, we shall compare the obtained poles with the poles of each individual 
publication, and draw certain conclusions about features of individual methods as well. 



C. Data base 



We start with a collection of data in which one part is fully available in the literature 
numeric values for the second part are provided by the authors (private communication refs. 

We have analyzed the following PWA amplitudes: 

1. Karlsruhe- Helsinki (KH80) [H - ttN elastic; 

As the influence of inelastic channels is in KH80 formalism introduced through forward and fixed cms 
scattering dispersion relations, KH80 docs not offer any inelastic channel amplitudes to be fitted. However, 
as we know that inelastic channels are extremely important in CMB formalism to ensure the stability of 
solutions (see following chapter and ref. H3), we have decided to constrain elastic KH80 amplitudes 
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with 7rN — > 77N WI08 amplitudes which fairly correctly depict the world agreement of the 77N channel 
amplitudes at lower energies - see Fig. [3] 

2. VPI/GWU - ttN clastic and ttN -> 77N; 

We have used single energy solutions (GWU-SES) [36| . and energy dependent solutions (WI08) [26l[36|. 

3. Giessen [3l[ - 7rN elastic and 7rN — > 77N 

4. EBAC - We have used two sets of PW solutions. Single-channel fit solution (7rN elastic fitted) — EBAC07 
[37| . and two channel fit solution (ttN clastic and ttN -» ?yN fitted) — EBAC08 with the ttN -> nN 
normalization adjusted in accordance with M. Doring and B. Diaz [4l|. with M. Doring and B. Diaz [4l[ 

5. Julich Hi| - ttN elastic and ttN -> 77N 

6. Dubna-Mainz-Taipei (DMT) [H, [H - ttN elastic and ttN -> rjN 

D. Fitting procedure 

We have used three channel CMB formalism with 7rN and 77N physical channels, and the third, effective 
two body channel to account for unitarity. We start with a minimal number of bare poles, and increase their 
number as long as the quality of the fit, measured by the lowest reduced x 2 value, could not be improved. In 
addition, a visual resemblance of the fitting curve to the data set in totality was used as a rule of thumb; i.e., 
we rejected those solutions that had a tendency to accommodate for the rapidly varying data points, regardless 
of the x 2 value. When the optimal number of poles is reached, we claim that we have found all partial wave 
pole solutions given by the chosen data set. As our criteria (minimal reduced \ 2 value and visual resemblance) 
are not extremely rigid, we have to differentiate between the two categories of poles: those which are seen with 
almost complete certainty, and those which are only consistent with the chosen set of data. The poles whose 
addition significantly improve the reduced % 2 value fall into the first category, those which improve the reduced 
X 2 value only marginally fall into the second one. It is interesting to note that in the latter case a number of 
almost equivalent, indistinguishable solutions for the questionable pole may be found. 

III. RESULTS AND DISCUSSION 

The intention of this article is to use only one method, Zagreb realization of CMB model, to extract 
pole positions from a world collection of partial wave data and partial wave amplitudes. As a test case, 
we do it for the Sn partial wave only. We use three channel model, with two measured channels 7rN, 77N, 
and the third channel 7r 2 N, which effectively represents all other inelastic channels, and "takes care of " unitarity. 

We extract pole positions from all available PWA and make a comprehensive analyses. We analyze the 
number of poles needed for a given partial wave, we discuss the importance of inelastic channels. 

A. Methodology 

The main feature of CMB multi-resonance, multi-channel model is a good control over determining the 
number of bare poles, and deducing the importance of number of fitted channels. 

Importance of inelastic channels 



The elastic 7rN scattering channel is the best measured and the most confident channel, so in all cases it 
is the pillar of the obtained partial wave amplitudes. Most of the information about the energy dependent 
structure of all solutions is coming from this channel, and it is expected that corrections are coming from other 
channels. Therefore, it is carrying the heaviest weight for obtaining final results. 

At this point we are bound to address one specific point in more details. 
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In ref. [4Cj we have discussed the continuum ambiguity problem in coupled-channel formalisms. Namely, 
once the inelastic channels are opened, it turns out that the differential cross sections themselves are not 
sufficient to determine the scattering amplitude. Let us illustrate why. If differential cross section da/dfl is 
given by |F| 2 , then the new function F = e i<s> F gives exactly the same cross section. It should be remarked 
that this phase uncertainty has nothing to do with the non-observable phase of wave functions in quantum 
mechanics. The asymptotic wave functions at large distances from the scattering center may be written as 
$f(x) w e l ' k ' x + F(6) s ^, r —> oo, so the phase of scattering amplitude is the relative phase of the incident and 
scattered wave. This phase has observable consequences in situations where multiple scattering occurs, and the 
continuum ambiguity is created. In the elastic region, the unitarity relates real and imaginary parts of each 
partial wave, and the consequence is a constraint which effectively removes this "continuum" ambiguity, and 
leaves potentially only a discreet one. The partial wave must lie on the unitary circle. However, as soon as the 
inelastic threshold opens, unitarity provides only an inequality: |1 + 2 i Fi\ 2 < 1 ImF; = |F;| 2 + //, where 

= t(1 — e~ 41mSl ). Therefore, each partial wave must lie upon or inside its unitary circle, and not on it. A 
whole family of functions $, of limited magnitude but of infinite variety of functional forms which satisfy the 
required conditions, does exist. However, in spite that they contain a continuum infinity of points, they are 
limited in extent. Thus, the islands of ambiguity are created. 

In ref. [4(| we have shown that including inelastic channels into the analysis is a natural way for eliminating 
continuum ambiguities. We have concluded that, by fitting only elastic channel, some of the resonant 
states which dominantly couple to inelastic channels might remain unrevealed, and we had to fit as many 
channels as possible. In the present paper we apply the following strategy: we shall first fit elastic chan- 
nel only, and show the poles we reveal. Then, we shall repeat the fit by fitting two channel processes, 7rN 
clastic and 7rN — > r/N data when available, and see how the number of poles, and their quantitative values change. 

The problem we are facing is the low quality input for the 77N channel, because 7rN — > 77N partial waves are 
in principle not well known. Anyway, as a final result, we have to accept the solution for which both channels 
are reasonably well fitted despite the low quality of the 77N channel data. 

Determining the optimal number of poles 



In CMB formalism the number of poles is a starting parameter. That in practice means that when fitting, we 
start with a minimal set of poles: one resonant and two for the background. Then we increase the number of 
resonant poles until the satisfactory fit is achieved, i.e., until the quality of the fit, measured by the reduced \ 2 
value, could not be improved. In addition, a visual resemblance of the fitting curve to the data set as a whole 
is used as a rule of thumb: we reject all those solutions which have a tendency to accommodate for the rapidly 
varying data points regardless of the x 2 value. 

In such a way we estimate the number of bare poles needed by our model, what in most cases corresponds to 
the number of resonant states. Observe that this is not so for dynamic resonances, i.e., for the dressed resonant 
states which do not have a corresponding bare pole. Therefore, what we compare is not the number of bare 
poles, but the number of dressed ones. (For a more extensive discussion on dynamic resonant states in 
Zagreb CMB model see ref. [3.) 

B. Fits 

We first fit 7rN elastic channel only. In accordance with the afore considerations, we first want to determine 
which resonances are we well determined only by this channel, and later on we want to see how much the 
inclusion of 77 N channel will modify the obtained result. 

1. tyN elastic channel only 
We show the result of the fit in Table HI The quality of the fit is shown in Fig. [5J 




FIG. 2: Agreement of 3R and 4R CMB curves with "world input" for single-channel fit. 
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TABLE I: "World collection" of poles for the single-channel fit, three and four resonant case. 



ANALYSES 



Fitted 
channel 



Number 
of 

resonances 



Bare poles 



W ai W S2 W S3 W S4 



MeV 



Dressed poles 



ReW 
-2ImW 



ReW \ / ReW 
-2ImW/ \-2ImW 
MeV 



ReW 
-2ImW 



Kb. 



KH80 



nN -)• nN 



1516 1638 1880 



1 188 1656 1713 2266 



1513 ' 

71 
1513 1 

113 



1661 
148 
' 1670 1 
194 



2263 
138 



WI08 



nN -)• nN 



1481 1657 3767 



1513 1624 1686 2517 



1492 1 

89 
1495' 

105 



1646' 

95 
1647 1 

81 



2396 
139 



GWU-SES 



nN -> nN 



1514 1645 2919 



1517 1650 1928 3768 



1500 ' 

106 
1505 ' 

97 



1646 1 
119 
1651 ' 
119 



2598 1 

210 
1944 1 

74 



2633 
345 



GIESSEN 



nN -> nN 



1464 1616 1731 



1474 1635 1718 2674 



1484' 

82 
1482 1 

82 



1641 
65 
' 1642 l 
65 



1861 
811 
1851' 
456 



2249 
287 



JUELICH 



nN -)• nN 



1518 1656 2177 



1511 1636 1719 2241 



1528 ' 

95 
1516 ' 

121 



1653 1 
110 
1654 1 
118 



2335 1 
372 
1665 1 
411 



2335 
403 



EBAC07 



nN -> nN 



1466 1641 2518 



1483 1643 1702 2237 



1498 1 
123 
1502 ' 
139 



1641 
89 , 

1638^ 
81 



2215' 

767 
1700 1 

408 



1862 
691 



EBAC08 



nN -> nN 



1515 1673 1826 



1512 1667 1980 3784 





2189 
637 



DMT 



nN -> nN 



1495 1643 2047 



1507 1647 1850 2100 



1486^ 

81 
1508 ' 

139 



2100 
212 



Averages 



1498(16) 

96(20) 
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114(20) 



1649(9) 
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1651(11) 
109(42) 



'2194(324) 
424(324) 

' 1793(108) ' 
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2253(224) 
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2. nN elastic and ttN — > rjN data 



As we have already mentioned, 7rN — > 77N data are rather old and vague, so the corresponding partial waves 
are poorly determined. Anyway, each analyzed PWA solution of our "world collection", with the exception 
of KH80 does offer some results for that channel, and we have consistently used it in the two channels fit. 
The only exception, KH80 amplitudes, do not have a corresponding 77N channel. We have been tempted to 
omit KH80 amplitudes from the coupled-channel analysis, but due to its extremely good analytical constraints, 
we have decided to keep it in some form. Instead of KH80 rjN channel, we have used the WI08 VPI/GWU 
solution believing that the Sn 77N channel amplitudes are confidently well known in the energy range s < 3 GeV 2 
(Tiab < 800 MeV), and in that range the WI08 VPI/GWU solution is a good numeric representation of a "world 
collection average" . 

We show the result of the fit in Table Ql] The quality of the fit is shown in Fig. [3] 



TABLE II: "World collection" of poles for the two channels fit, three and four resonant case. 



ANALYSIS 
(Fitted channels) 



Number 
of 

resonances 



Bare poles 



W S1 W S2 W S3 W S4 
MeV 



Dressed poles 
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-2ImW 
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MeV 
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_> 
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1517 1637 1865 
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1670 
163 
1650 ' 
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Averages 
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1986(373) 
471(301) 
1745(80) 
220(95) 
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FIG. 3: Agreement of 3R and 4R CMB curves with "world input" for two-channel fit. 
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All obtained pole positions are shown in Fig. [4] 
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FIG. 4: Poles of a "world collection" of PWA. 
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C. Individual comparison 

Preliminary considerations 

As it has been generally accepted, T-matrix pole positions are the most recommendable singularities to 
be compared with QCD. However, obtaining them definitely means going into the complex energy plane while 
having at ones disposal only the physical T-matrix values (values for the real energy) . This analytic continuation, 
however, has to be a model dependent procedure by definition, because there is no a priori rule how to choose 
the analytic functional form which is to represent a measurable subset out of all possible T-matrix values. 
Therefore, the reader has to be fully aware that the pole positions we find, and the pole positions given by 
the original publications have to be different by definition, and the reason is that each investigated "world 
collection" solution has its own way how to analytically continue the measurable physical T-matrix values. 
However, comparing the number of needed poles, their distribution and genesis (genuine or dynamic) obtained 
by our approach with those from original publication is certainly justified. It is also a convenient way to 
establish whether a certain pole is only a result of a poor knowledge of measured process, or indeed is a genuine 
singularity needed by the data, but still not yet well established. So, hereafter, we analyze qualitative features 
of the partial wave singularity structure, and intentionally avoid to compare their numeric values. 

1. KH80 

The KH80 amplitudes are essentially single-channel partial wave data with some information about inelastic 
channels introduced through forward dispersion relations, and analyticity strictly imposed on the level of 
fitting procedure using Pietarinen expansion [42j. As no assumption on the analytic functional form about 
partial-wave amplitudes has been done, search for resonance parameters is a separately defined procedure. 
Brcit-Wigner parameters are obtained as a local fit in the resonance region with background contribution 
unitary added on the level of S matrices, and poles are extracted using single-channel pole position extraction 
methods (speed plot and Argand diagram). Original publication reported two poles. 

In our approach we concur the existence of first two poles, and we find them strongly dominated by the 
elastic channel. 

The third N(2090) pole is in our fit definitely needed. The fourth pole is allowed by our fit in both, single 
and coupled-channel constellation (improvement of the reduced xfj), but its quantitative constraint will need 
more inelastic channels than only 17N. In each configuration numerical values of third and fourth pole are not 
yet sufficiently well constrained. 

2. GWU-SES and WI08 

While the original publication gets the pole positions by analytically continuing energy dependent solution 
into the complex energy plane, an obvious advantage of our approach is that we can obtain the pole positions 
independently from both, single energy (GWU-SES) and energy dependent (WI08) VPI/GWU solutions. We 
have to remember that VPI/GWU pole positions are extracted from the analytic form determined by their 
Chew-Mandelstam K-matrix approach, which is fitted directly to the data, and not to their single channel 
solutions. Consequently, the pole positions "corresponding" to their single energy solutions are by them not 
yet discussed. In this paper we may use the same formalism for both, single energy and energy dependent 
solutions, and treat them as an independent input. Hence, we get two sets of solutions. 

The general conclusion for both VPI/GWU solutions is the same, and it is very similar to the findings for 
the KH80 input: we confirm the existence of first two poles, and find them strongly dominated by the elastic 
channel. The third N(2090) pole is in our fit definitely needed. The fourth pole is allowed by our fit in both, 
single and two channels constellation (improvement of the reduced Xr)j but its quantitative constraint will 
need more inelastic channels that 77N. 
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It is very interesting to compare WI08 with GWU-SES. In spite of the fact that the WI08 solution is seemingly 
very smooth above the second peak, definitely much smoother than the GWU-SES solution, our model still 
requires the third and fourth pole almost in a same way for both solutions. The need for a third and fourth pole 
for the smooth WI08 solution came as a surprise for us. Quantitatively, all pole positions are similar for both 
solutions: quite well defined for the first two poles, dominantly determined with the elastic channel. Inclusion of 
inelastic rjN channel data modifies first two pole positions only slightly. However third and fourth pole positions 
remain strongly influenced. 



3. DMT amplitudes 



DMT collaboration has originally looked for the pole positions using the speed plot technique. They 
have established the existence of three poles, N(1535) , N(1650) and a third pole corresponding to N(2090) 
- see ref. (30| . However, triggered by their old research of photo-production channels in which they had 
established the strong probability for the existence of new S-wave resonant state in the vicinity of 1846 MeV 
29, 32|, they have recently repeated the analysis and confirmed the existence of this new state at 1880 MeV 17 1. 



It is interesting to note that our procedure for DMT amplitudes also indicates the existence of 4 poles. As 
seen in Fig. [2J our 3-resonant fits do miss some structure in elastic partial waves at higher energies requiring 
the increase in the number of parameters. Repeated fits with 4 resonances rectify this problem and at the same 
time show a significant improvement in the reduced % 2 . So, our fits concur their latest findings that the DMT 
Sn solution really contains 4 poles 



4- EBAC amplitudes 



EBAC has produced three sets of partial wave amplitudes: the first, single channel set where only 7rN elastic 
data have been fitted - EBAC07 - [331, an d two additional sets of amplitudes where data from more than one 
channel was used to constrain the fit; in t his p articular case 7rN and 77 N channels. The unpublished set [43j in 
a way supersedes the former 2008 analysis |38| where unpleasantly large change of 7rN elastic partial waves was 
needed to accommodate for the second channel. We have analyzed both sets of amplitudes wandering whether 
a significant change in poles between the two is found. However, as no numeric data for the unpublished set 
is available to us, we have attempted to "read off" the data directly from the graph, and that has introduced 
uncontrollable numeric instabilities. Therefore, we have decided to omit the EBAC10 preliminary data from 
our analysis until the final results are published. 



The EBAC group has in all three analyses used two bare poles, situated relatively high in energy (M > 1.8 
GeV), and reported two dressed poles corresponding roughly to N(1535) and N(1650). Third and fourth pole 
have not been found. Just as a preview, we can state that our analysis finds all three solutions very similar. For 
all three sets we confirm the existence of the first two poles, and they are strongly constrained by the 7rN elastic 
channel alone. However, our fits indicate that significant improvement reduced x 2 is achieved if the third and 
fourth poles are allowed for. These poles are needed by the fit, but still poorly determined by only two inelastic 
channels. 



5. Jiilich amplitudes 

Similarly to many, Jiilich group fits their model to VPI/GWU data (to energy dependent WI08 set (HI), 
and very much like WI08, obtains a very smooth behavior above 1800 MeV. The only difference with respect 
to WI08 is a different behavior of high energy tail: while the real part of Jiilich amplitudes falls with energy 
and the imaginary part raises, in case of WI08 amplitudes the result is just opposite. Therefore, a difference 
between the two should not be found in cross section measurements, but only possibly in some polarization 
ones. They also report two Sn poles. 

Consequently, we expect that our results for pole positions of Jiilich amplitudes show a very similar behavior 
to WI08, and that is fulfilled. 
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The most prominent feature of our analysis of WI08 amplitudes — that in spite of smooth high energy behavior 
we need more than two poles to fit the input — is confirmed for Jiilich amplitudes as well. It is completely clear 
that we need at least three poles to satisfactorily reproduce the amplitude shape, and their amplitudes are 
in our analysis consistent with four Sn poles. Very similar as before, the third and fourth poles are rather 
undetermined with only two channel constraints. We have discussed the possibility of finding extra poles in 
Jiilich amplitudes with M. Doring in Zagreb last fall fill ], and this possibility has not been entirely ruled out 
even in analytical continuation Jiilich method. They have simply not looked for the pole in that energy range. 
However, even while this pole might be around 1800 MeV, it must be rather far in the complex energy plane. 



Giessen group also fits GWU-SES data in 7rN elastic channel, and gets a reasonable agreement with the 
input. The main difference with respect to "world collection" amplitudes, again lies in the 77N channel data. 
Nost results for this channel more or less agree within the N(1535) dominance range, but significantly deviate 
in the higher energy region. 

The Giessen model assumes K-matrix Born approximation where the real part of the Green function is 
neglected, the analyticity is manifestly violated. Consequently, the comparison of poles obtained in our 
fit with poles of these amplitudes is more questionable, as the main assumption for the correct analytic 
continuation — that is the analyticity of the model — is not preserved for both models. 



As the main aim of the paper is to use one method in order to eliminate systematic uncertainties in pole 
extraction, we summarize our primary results. 



All pole positions and their averages are shown in Figs. [5]and|ni 

Three resonant case 

As the number of accepted Sn resonances in PDG 19] is three, we first stooped our fit at 3 bare 
poles. 

By inspecting 3R solutions in Table 1 and Fig. [5] we observe: 

• First two poles N(1535) and N(1650) are extremely well determined in all PWA. 

• We find their average value to be 



• All PWA do need a third pole, but its position is extremely ill-defined; KH80, Giessen and EBAC08 prefer 
the values between 1700 and 2000 MeV, while the rest have the values above 2000 MeV. 

• The resulting average value is poor 



6. Giessen amplitudes 



D. Primary result: Averages 
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FIG. 5: "World collection" of poles for the three resonance single-channel fit. 



This separation in two preferred ranges of the third pole among different PWA permits us to speculate whether 
the fitting rules allow for the existence of the 4-th pole. 

Four resonant case 

We have repeated the fit with 4-bare poles, and results are collected in Table 1. as 4R solutions. We 
show the result in Fig. [6] 
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FIG. 6: "World collection" of poles for the four resonance single-channel fit. 



By inspecting 4R solutions in Table 1 and Fig. [5] we observe: 



• We have found that all other PWA if not require, then are at least consistent with the 4 Sn poles, even 
the EBAC amplitudes which are based on only 2 bare poles input. 
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• The reduced \ 2 is either improved, or at least stays the same for all solutions; that justifies the inclusion 
of the fourth pole. 

• First two poles N(1535) and N(1650) are again very well determined in all PWA 

• We find their average value to be 



• Contrary to our expectations, and in spite of the fact that the reduced \ 2 is improved practically every- 
where, the scatter in 3rd and 4th pole remain. 

• The resulting average value for third and fourth pole is poor 



The existence of the fourth pole is not convincing. 

Due to the fact that third and fourth pole poorly couple to the elastic channel that is only used at this instant, 
we conclude that fitting other channels is inevitable if the improvement on the third and fourth pole parameters 
is to be achieved. 



The poor determination of third and fourth pole for the single channel fit confirms our former findings that 
inelastic channels are essential for fully constraining all resonant states (scattering matrix poles) - see ref. 
poj . The problem with stability of minimization solutions lies in the fact that the 77N channel data are old, 
scarce, and unreliable (for instance Brown data at higher energies - see discussion in ref. [HI), so 77N channel 
partial waves are imprecise. Even when being of lower quality, the rjN channel data still represent a valuable 
constraining condition, because the general trends of 77N channel are to be simultaneously reproduced together 
with the details of clastic channel, and that is by no means simple. The results of the fit are given in Table 2 
and Figs. [7] and H 

Three resonant case 

As the number of accepted Sn resonances in PDG 1S| is three, we first stooped our fit at 3 bare 
poles. We show the result in Fig. [7] 

By inspecting 3R solutions in Table 2 and Fig. [7] we observe: 

• First two poles N(1535) and N(1650) are extremely well determined in all PWA. 

• We find their average value to be 





2. tyN elastic and -rN — > nN data 




• If we compare these numbers with the result of single-channel, three resonance fit: 
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FIG. 7: "World collection" of poles for the three resonance, two channels fit. 



M1535) Sii = ^ 5S±i6 

Ar(l 650 ) 5ll = f^ 5±9 
\J0J ± 25 

we see that the difference is within one standard deviation. Real parts of the resonances are 
almost completely reproduced, while the imaginary parts are slightly shifted downwards. 

• All PWA do need a third pole, but it's position is again extremely ill-defined. 

• The resulting average value is poor 



A(2090) Su 
Four resonant case 



1986 ± 373 
471 ± 301 



We have repeated the fit with 4-bare poles, and results are collected in Table 2. as 4R solutions. We 
show the result in Fig. |8] 

By inspecting 4R solutions in Table 1 and Fig. [5] we observe: 

• We have found that all PWA if not require, then are at least consistent with the 4 Su poles, even the 
EBAC amplitudes which are based on only 2 bare poles. 

• The reduced x 2 is either improved, or stays the same for all solutions. That justifies the inclusion of the 
fourth pole. 

• First two poles N(1535) and N(1650) are again extremely well determined in all PWA 

• We find their average value to be 
'l502 ± 23^ 



JV(1535) Sn = 
A(1650) Su = 



144 ± 39 

1649 ± 13 
117 ±37 
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FIG. 8: "World collection" of poles for the four resonance, two channels fit. 



The resulting average value for third and fourth pole are: 
N(xxxx) S n = ( 1745±8( A 

iV(2090) S n = f2i9i±24iy 
v ' \ 392 ± 301 J 

The scatter in the 3rd pole is significantly reduced, and the indications for its existence are strong. 

The existence of the fourth pole is strongly indicated, but still not quite convincing. 

If we compare these numbers with the result of single-channel, four resonance fit: 

JV(1535) Sn = (^02±12\ 
\ 114 ±20 J 

7V(1650) Sn = f 1651 ^ 11 ), 
\ 109 ±42 J 

1 r F7 r-p— (l793±108\ 

\ 453 ± 325 J 

7V(2090) Su = ( 2253 ±22 A 
\ 356 ± 212 J 

we conclude: 

the 77N channel data have confirmed the good constraint on N(1535) and N(1650) Sn states, they have 
improved the confidence limits for the existence of the new N(xxxx) Sn state, but they are definitely 
insufficient to constrain the fourth Sn pole. 

Therefore, other channel partial waves have to be included. 
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IV. CONCLUSIONS 



We have offered one model, the Zagreb realization of CMB model, for extracting pole positions from 
a "world collection" of partial-wave amplitudes which we treat as partial wave input data, and extracted 
the results. Using only one method enables us to make a statistical analysis of partial wave poles in a 
manner that we avoid the systematic error caused by the different assumptions on the partial wave ana- 
lytic function form. We have in details explained the idea, and presented the results for the Sn partial wave only. 

We have analyzed the single channel fit (only one channel data are used to constrain the fit), and in details 
investigated what are the consequences of enlarging it to a two-channel one with the 77N channel. We have 
concluded that even low quality data in the second channel are sufficient to notably constrain the arbitrariness 
of the poorly determined poles. However, we also concluded that for the third and fourth S-wave poles, 77N 
channel is not sufficient. 

We found that the first two Sn poles are extremely well defined by elastic channel and that the included 
inelastic 77 N channel introduces only small modifications of the elastic channel result. 

We have shown that all members of the partial wave "world collection", in spite of the fact that some of 
them have assumed only two S-wave resonant states, are consistent with at least three T-matrix poles. We 
have also demonstrated that there is a strong statistical indications that the fourth pole is present in each of 
the "world collection" member, despite the fact that no one has seen it up to now. The exception is the latest 
DMT analyses [l7| where all four poles are detected. The new Sn pole we have found can be identified with 
the Sn(1846) pole seen in photo-production channel by the DMT collaboration j29l |32|. 

We affirm that the results of the 4-resonant, double channel fit should be treated as a final result, and we 
offer the world average: 




(12) 
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